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ABSTRACT 

Wc follow the growth of baryonic structure in the presence of a magnetic field within 
an approximate cosmological magneto-hydrodynamic simulation, produced by adding 
an (isotropic) magnetic pressure related to the local gas pressure. We perform an en- 
semble of these simulations to follow the amplification of the field with time. By using 
a variety of initial field strengths and changing the slope of the power law that governs 
the way the field grows with increasing density we span the range of current observa- 
tions and demonstrate the size of the effect realistic magnetic fields could have on the 
central density of groups and clusters. A strong magnetic field significantly reduces 
the central gas density which, in turn, reduces observable quantities such as the X-ray 
luminosity. 

Key words: magnetic fields, hydrodynamics, methods: numerical, galaxies: clusters: 
general 



1 INTRODUCTION 

' In an era where we have such a firm grasp of many 
of the fundamental parameters that define the universe 

■ in which we live (Spergel et al. 2003), the underlying 
framework of cosmology based on gravitational insta- 
bility of cold dark matter is no longer subject to seri- 
ous doubt. However, although the idea that large po- 
tential wells are built up by the hierarchical assem- 

■ bly of many thousands of smaller sub-blocks (Press & 
Schechter 1974; Davis et al. 1985) has not been seriously 
challenged for decades, there is much current debate 
about the presence or absence of large amounts of sub- 
structure (Kazantzidis et al. 2004) and about the cuspi- 
ness of the inner parts of dark matter profiles (Moore 
et al. 1999; Power et al. 2003; Diemand, Moore & Stadel 
2004). Moreover, the reaction of baryonic material to 
this underlying framework is not well understood at all. 
We still do not readily understand why such a small 
fraction of baryons cools to form stars (Fukugita, Hogan 
& Peebles 1998; Balogh et al. 2001), how the first stars 
are distributed (Abel, Bryan & Norman 2002) and sub- 
sequently reionize the universe (Sokasian et al. 2004) or 
why the gas at the centre of dark matter haloes obsti- 
nately refuses to cool efficiently (Ponman et al. 1999; 
Pearce et al. 2000, 2001; Voit & Bryan 2001). 
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An additional complication arises from the well- 
known fact that the intra-cluster medium (ICM) in 
galaxy clusters is magnetized. This is inferred from 
observations of diffuse radio haloes and hard X-ray 
emission as well as Faraday rotation measurements 
(Kronberg 1994; Feretti & Giovannini 1996; Govoni et 
al. 2001a; Taylor et al. 2001; Eilek & Owen 2002; Murgia 
et al. 2004; Fusco-Femiano et al. 2004). Magnetic fields 
are observed throughout the Universe at a wide range 
of epochs and scales; from a few \iG up to a few mG 
within galaxies (Beck 2000; Krause 2003), and from a 
few tenths of a /xG on cluster and super-cluster scales, 
to values in the centre of cooling flow clusters reaching 
as high as a few tens of \iG (Carlstrom, Holder & Reese 
2002). At high redshifts magnetic fields of a few nG have 
been observed in systems such as Lyman-a absorption 
systems, (Kronberg & Perry 1982; Norman 1990; Kron- 
berg 1994; Oren & Wolfe 1995), radio-galaxies (Athreya 
et al. 1998; Pentericci et al. 2000), and proto-clusters 
(Pentericci et al. 2000; Bagchi et al. 2002). Although 
magnetic fields of the order of /zG in strength may 
seem insignificant, the vast scales involved mean that 
a significant amount of energy is stored in these fields, 
and consequently they can be dynamically important. 
Where they come from is still a mystery. In this paper 
we shall assume that some primordial seed field exists 
and experimentally determine whether it can have a 
significant effect on cluster properties. 

Analytical and numerical studies provide an impor- 
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tant counterpoint to these observations, aiding our un- 
derstanding of the distribution and perhaps origin of the 
magnetic field. The first attempt to implement MHD 
(magneto-hydrodynamics) into SPH (smoothed particle 
hydrodynamics) was the polytropes studied by Gingold 
& Monaghan (1977). Further aspects of the applica- 
tion of SPH to magnetic phenomena were considered by 
Phillips & Monaghan (1985), Phillips (1986) and Mon- 
aghan (1992). Phillips & Monaghan demonstrated that 
when the MHD equations are written in a conservative 
form a numerical effect occurs: SPH particles tend to 
clump due to the presence of an artificial tension. Sub- 
sequent authors (Meglicki, Wickramasinghc & Dewar 
1995; Cerqueira & de Gouveia Dal Pino 2001; Hosking 
2002) circumvented this problem using non-conservative 
forces but this approach does not work well when shocks 
are present. More recently, Price & Monaghan (2004) 
developed a new approach in which they added a small 
artificial stress that prevented the numerical instability. 

Detailed MHD simulations of the growth of struc- 
ture in the universe pointed out the importance of 
merger events on the distribution and strength of the 
final magnetic field and that the initial field structure is 
completely wiped out during cluster formation. In addi- 
tion, compression and shear flows can strongly amplify 
the seed field (Birk, Wiechen & Otto 1999; Roettiger, 
Stone & Burns 1999; Dolag, Bartelmann & Lesch 1999, 
2002). 

Loeb & Mao (1994) have suggested that, if the re- 
sultant magnetic field were sufficiently strong, the ICM 
could be supported to a significant extent by magnetic 
pressure in addition to thermal gas pressure. In clus- 
ters where this was the case the magnetic field would 
be dynamically important (Dolag et aZ. 2001; Eilek & 
Owen 2002), so this support might contribute to the 
discrepancy between X-ray and gravitational lensing 
mass estimates of the central regions of galaxy clusters. 
Goncalves & Friaca (1999) found that the magnetic field 
could be dynamically important on scales as small as 
< 1 kpc, although it seems unlikely that either these 
small scale fields or a more widely distributed field such 
as that proposed by Loeb & Mao could be the main 
reason for the discrepancy between the mass estimates 
in the central regions, at least in relaxed clusters (Dolag 
et al. 1999; Dolag & Schindler 2000). 

In this paper we examine the effect of introduc- 
ing a large-scale magnetic field on the growth of bary- 
onic structure. As gas collapses isotropically the mag- 
netic field strength grows as B oc assuming that 
the field lines are frozen into the plasma. Accordingly, 
since the magnetic pressure is P mag oc B 2 , the addi- 
tional gas pressure due to the presence of a magnetic 
field rises as p~s . Here we implement this simple assump- 
tion into an N— body hydrodynamics code and follow 
what happens to the density of each of the gas parti- 
cles. Obviously this method relies on a smooth collapse 
of the gas without massive tangling of the magnetic field 
lines, a situation that restricts its applicability to rela- 
tively low values of the overdensity. For this reason we 
have made no attempt to follow gas cooling or galaxy 
formation, but have concentrated on the hot haloes of 



galaxy groups and clusters. In order to mimic more real- 
istic fields we parametrise our model so that the power- 
law index for the B — p relation is one of the free pa- 
rameters, a; the motivation for this is discussed in the 
next section, where we detail our numerical models and 
magnetic field approximation. Our results arc presented 
in Sections 3 & 4, where we study various initial field 
strengths and power law dependencies with density as 
well as varying the redshift at which the field is initiated. 
Section 5 contains a discussion and our conclusions. 



2 NUMERICAL CALCULATIONS 
2.1 The Simulation Method 

The basis for our model of the passive evolution of mag- 
netic fields is Hydra, an AP 3 M - SPH code (Couch- 
man, Thomas & Pearce 1995). Smoothed-particle Hy- 
drodynamics (SPH) is a Lagrangian numerical method 
which follows the motion of a set of fluid (gas) elements 
represented by discrete particles. The thermal energy 
and velocity of each particle are known at any given 
time and each particle has a fixed mass. Properties of 
the gas at the position of a particle can be estimated 
by smoothing these quantities over the -/Vsph nearest 
neighbouring particles. The gas properties are then used 
to calculate the forces acting on each particle in order 
to update the positions and velocities. In cosmologi- 
cal simulations both dark matter and gas particles are 
included and the particles are initially distributed in a 
manner consistent with a cosmological power spectrum. 
If the process of galaxy formation is to be simulated 
then radiative cooling of the gas must also be included; 
we neglect this for the purposes of this paper. 

The base simulation used throughout this paper has 
128 3 gas and 128 3 dark matter particles with individ- 
ual masses of 6.58 x 1O 8 /i" 1 M and 4.27 x 10 9 /i _1 M Q 
respectively, in a periodic box of 50/i _1 Mpc. The power 
spectrum is that appropriate to a cold dark matter uni- 
verse with the following parameter values: mean mass 
density parameter Qm — 0.3, cosmological constant 
Oa = 0.7, baryon density parameter Slf, = 0.04, Hub- 
ble constant (in units of 100kms _1 Mpc _1 ), h = 0.7, 
power spectrum shape parameter T — 0.21 and rms 
linear fluctuation amplitude er 8 = 0.90. These cosmo- 
logical parameters are close to those obtained through 
fits to WMAP data (Spergel et aZ.2003): fl M = 0.27, 
A = 0.73, n h = 0.044, h = 0.71, a s = 0.84. The gravi- 
tational softening length is 20/i~ 1 kpc, fixed in physical 
coordinates. The starting redshift was z s tart = 49. 

We carried out a set of 9 runs, each one charac- 
terised by different parameter values (see Table 1 and 
following section). For each run we extracted clusters 
with a spherical overdensity criteria and selected those 
with at least N = 1000 particles of each kind within 
the virial radius. The number of objects we analysed is 
of the order of 140 per run. Given the limited volume 
of our simulations, we have only a couple of moderately 
rich objects, with M vir > 10 14 /i _1 M Q . 
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Figure 1. Magnetic field strength as a function of overdensity. 
The shaded regions and the continuous line delimit our range of 
a and Bq. Top green: Bo = 10~ 9 G and 0.5 < a < 1.; bottom 
blue: Bo = 10 — 10 G and 2/3 < a < 1.1; red continuous line: 
So = 10 — n G a = 0.8. The black dashed line is from Dolag et al. 
2005. 
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Table 1. Parameters for the simulations used in this paper. By 
column a simulation identifier, the magnetic field normalization in 
Gauss (-Bo) an d the power law index of the field— density relation 
(a). 

2.2 Approximating a Magnetic field 

In this section we describe how we introduce a mag- 
netic field into our simulations and justify our choice of 
parameters. 

According to the standard picture of magnetic field 
evolution a seed field generated at high redshift is am- 
plified by adiabatic compression and merger events oc- 
curring during cluster formation. A straightforward an- 
alytical argument shows that under the assumptions of: 

• B-field frozen into plasma 



• uniform spherical collapse 

• small magnetic field 

• mass conservation 

• magnetic flux conservation, 

then B scales with the gas density p as B oc p a with 
a = 2/3. On the other hand, if the magnetic field is 
important, there will be a preferential axis x for the 
collapse due to the presence of strong field lines and 
the collapse proceeds cylindrically. This, together with 
flux conservation and the assumption that gravity is in 
equilibrium with thermal pressure along the symmetry 
axis (2nGpx 2 ~ c 2 s ) leads toBoc (Tp)? , so if the object 
is isothermal the power-law relation between B and p 
acquires a different value of a = 1/2 (Crutcher 1999). 

More detailed numerical studies have found that in 
realistic large-scale structures without special geome- 
tries the amplification of the seed field is not quite as ex- 
pected from these simple collapse models, with a rough 
power-law dependence but a higher value of a (Roet- 
tiger et al. 1999; Dolag et al. 1999; Dolag et al. 2005). 
King & Coles (2006) showed that an average over small- 
scale collapses leads to an average value of a which is 
higher than the isotropic value a — 2/3; the average 
of a non-linear function is not the same as the non- 
linear function of the average. Moreover, other phenom- 
ena such as mergers and anisotropics act to enhance 
the field still further. It is therefore likely that different 
power indices better describe the field behaviour in dif- 
ferent density ranges, as found observationally: Dolag 
et al. (2001) found a = 0.9 for the galaxy cluster Abell 
119, while Crutcher (1999) got a = 0.47 for molecular 
clouds in a density regime of about 10 3 — 10 4 cm~ 3 . 

The parametrisation we have chosen corresponds to 
a polytropic equation of state of the form P cx p 2a ; large 
values of a will therefore produce a steeper relation than 
pertains for the thermal pressure (a ~ 5/6). This is the 
regime in which interesting physics can occur. 

These considerations all lead us to approximate the 
magnetic field B(z) as: 

B(z) = B (l + z) 3a p a , (1) 

where p is the local overdensity, Bq the field in an un- 
perturbed universe, both in comoving units, and z is the 
redshift. The term p a accounts for the field amplifica- 
tion by incorporating a as a free parameter. By varying 
the power-law index we emulate the different possible 
amplification mechanisms, though we stress that we do 
not solve the full MHD equations exactly, in particular 
we ignore the back-reaction of the field on the gas which 
would break the simple scaling of equation (1). In real- 
ity it is likely that the amplification of a -B-field by the 
factors we have discussed would saturate around the 
cquipartition value. Care must be taken, therefore, in 
interpreting our results when magnetic pressure domi- 
nates the thermal pressure. We are primarily interested, 
however, in the growth of structures through the linear 
and quasi-linear regime so we hope our calculations have 
a reasonable regime of validity. It is also important to 
understand that we can not use this method to simu- 
late ordered fields displaying large-scale anisotropy: we 
are restricted to fields that are tangled on a sufficiently 
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small scale that their effect on the large-scale gas mo- 
tion is isotropic. 

The standard SPH momentum equation for particle 
i using a smoothing kernel W is: 



§ =-E^(§+f+n]v^, 



(2) 



and the corresponding equation for the rate of change 
of the internal energy e, is 



(3) 



§-;i>(M*n)« 

where z/y is the relative momentum of particles i and 
j. In these equations Pi is the pressure measured at the 
position of the i-th particle. 

Given the presence of a magnetic held, the magnetic 
contribution to the pressure can be written as 

Fmag = w (4) 

with /io = 4ir x 10~ 7 and £? is interpreted as the rms 
value of the (tangled) held. If we assume that the mag- 
netic held is passively moved and squeezed by the gas 
we can think of this magnetic pressure as being related 
to the baryonic density in the same way as the thermal 
pressure but with a different equation of state described 
by the parameter a. Writing P = P t h + -Pmag where Pth 
is the thermal pressure (calculated using the standard 
SPH methods) we can accommodate both forms in the 
standard SPH equations by simply adding an extra ef- 
fective pressure variable. Note, however that this ap- 
proach has severe limitations: there is no back-reaction 
of the held on the gas; it cannot cope with large-scale 
magnetic structures which are inherently anisotropic; 
and the polytropic model we use does not rigorously 
conserve energy. One should not place too much literal 
emphasis on the quantitative results we present, espe- 
cially for high densities. 

The magnetic held strength as a function of over- 
density at z = in our models is shown in Figure 1. 
We hxed the normalization factor Bq and the power 
law index a so that we span the range of observed held 
strengths in halo cores as well as the values obtained by 
the more complex model of Dolag et al. (2005) which is 
shown as the dashed line. The top green shaded region 
represents the area in the B—p plane where Bo — 10 _9 G 
and 0.5 < a < 1.; the lower blue region is relative to 
B Q = 10~ 10 G and 2./3. < a < 1.1 while the red contin- 
uous line is Bq — 10~ n G and a = 0.8. 

For these values the highest normalisation models 
challenge the latest observational constraints for the 
magnetic energy density in clusters obtained from Fara- 
day rotaton measures (Vogt & Ensslin, 2006). Our mod- 
els were chosen to span the full range of possible mag- 
netic held strengths, with the strongest helds larger 
than observed (and so producing consequences larger 
than expected in the observed Universe). 

In uncollapsed regions the observed strength of the 
magnetic held is currently only an upper limit. Sev- 
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Figure 2. Mean radial fractional pressure. Each line represents a 
different run with the parameters labelled and displays the mean 
ratio of the thermal to magnetic pressure in concentric spherical 
bins over all clusters. 
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Figure 3. A comparison at different redshifts of the densities of 
equivalent gas particles from simulations without a magnetic field 
(noB) and one with a magnetic field of strength 10 -9 G, a = 1. 
(B0960, left panel) and a magnetic field of strength 10 — 10 G, a = 1 
(B1060, right panel). The appropriate redshift is labelled at the 
end of each profile. 



eral of our power-law models exceed this value but as 
no structures collapse and we are not analysing these 
volumes this has no affect on our results. Our weakest 
held strength was chosen so as to fall within the range 
of allowed values in the voids. 
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Once the amplitude and power-law index of the mag- 
netic field have been chosen there remains the choice of 
the redshift at which the field is first imposed. If the 
magnetic field is turned on at late times, after much of 
the structure has formed, then it has little affect on the 
resultant objects. If, however, the field is imposed before 
significant structures form (which in turn depends on 
the resolution of the particular simulation being stud- 
ied) then the initial redshift makes little difference to 
the final distribution of the matter. 

Using trial simulations we established that for rea- 
sonable choices of initial field, z > 5 is necessary to 
have significant effect, but at higher redshifts the spe- 
cific choice of z is not important. We therefore impose 
the magnetic field at z = 9 for all subsequent runs. 



3 REDSHIFT EVOLUTION 

We start our analysis with Figure 2 where we show 
the ratio of P ma g versus Pth as a function of radii. The 
values have been obtained averaging over all clusters 
and each line is relative to a different run. The inclu- 
sion of a magnetic field B adds a non-thermal pressure 
Pmag oc B 2 that counteracts gravity. When P mag is not 
negligible compared to Pth the magnetic fields become 
dynamically important in that the additional magnetic 
pressure supports the gas against gravity and so for the 
same temperature the gas can reside at a lower density. 
Also, the additional pressure can slow the infall of the 
gas, retarding the evolution of structure. In reality we 
found that P mag is orders of magnitude smaller than 
Pth in most cases and the ratio of gas magnetic to ther- 
mal pressure generally falls with radius. However, even 
an additional 1% pressure support in the core region 
of collapsed objects can have serious consequences, par- 
ticularly for observables such as the X-ray luminosity 
which depends on the square of the gas density. 

In Figure 3 we plot the rank-ordered density of all 
the particles in the run with magnetic field compared 
to the rank ordered density without a magnetic field at 
a range of different redshifts. With the larger field (left 
panel) the high density particles incur an appreciable 
drop in their density. At lower field values (right panel) 
the effect of the magnetic field is much reduced. 

As can be seen for the high field redshift zero case, 
at low redshifts and high field values particles near the 
mean density actually end up denser than when no mag- 
netic field was present. This is due to the magnetic field 
pressure slowing the infall and reducing the strength 
of the accretion shock, resulting in a lower final entropy 
and consequently greater final density. We find that 70% 
of the particles with density ratio po960 / PnoB > 1 have 
entropy ratio so96o/snoB < 1, while the percentage in- 
creases to more than 90% if we consider particles with 
Po96o/PnoB > 10 (where p 096 o (so96o) and p noB (s no B) 
are the densities (entropies) for run 0960 and noB). 

In Figure 4 we show a series of gas density maps 
(where brighter colours indicates denser gas) of a box 
of side 10/i _1 Mpc centered on one of the most massive 
clusters (1.5 x IO^/i^Mq). In this sequence of 6 pan- 
els the field strength generally increases down the page 




Figure 5. Map of the ratio of magnetic to thermal pressure at 
z = for B1060. The box is 10h~ 1 Mpc across, same region as in 
figure 4. 



with the precise configuration of the magnetic field for- 
mulation that corresponds to each panel given in Ta- 
ble 1. The presence of a strong magnetic field also re- 
sults in a smoother mass distribution, with small clumps 
washed out and less frequent shocks. 

Spatially mapping the ratio of the magnetic to ther- 
mal pressure as in Figure 5 (which shows the same 
region as in Figure 4) reveals that only in small ob- 
jects and cluster cores can the magnetic pressure ap- 
proach equipartition with the thermal pressure (Pmag ~ 
0.1 — 0.3Pth), provided that Bq and a are sufficiently 
high (Bq > 10~ 10 and a > L). This results in a shal- 
lower density profile (Figure 6) and a lower entropy 
level in the inner regions of collapsed objects, while the 
outskirts are not affected. The mean radial gas density 
profile is displayed in Figure 6 which shows that field 
strengths not much higher than those observed can sig- 
nificantly depress the central gas density. Moreover we 
expect small structures to be affected more than big 
clusters because the former are characterised by a shal- 
low gravitational potential and so a relatively small ad- 
ditional pressure can be sufficient to deter gas infall. 
This can be easily seen from density maps (Figure 4) 
but, due to the small simulation volume, we do not have 
a wide range in mass to properly address any trend 
in mass from density profiles. Also, unfortunately, our 
spatial resolution is not high enough to properly in- 
vestigate the regions where the magnetic field effects 
should be more pronounced. In Figure 7 we show the 
bolometric luminosity in erg s^ 1 versus the bolomet- 
ric emission weighted temperature in KeV for differ- 
ent runs. The reduction in the central gas density low- 
ers the bolometric X-ray luminosity of a group-size ob- 
ject (M < lO^h^M©) by 40 - 70% for field values of 
B Q = IQ- W G and a = 1.1 (Figure 7). The reduction in 
luminosity can be even more than one order of magni- 
tude for the strongest field we analysed. 
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Figure 7. The luminosity— temperature relation. The values plot- 
ted are the emission weighted temperature and the bolometric lu- 

Figure 6. Density profile at z = for our clusters. The lines minosity within r 500 . Observations are from Ponman et at. 1996; 

represent the first and the third quartile, the dashed blue area is Helsdon et al. 2000 and Novicki et al. 2002. 

for run B1066 and the white dotted one for the noB run. 



4 OBSERVATIONAL CONSEQUENCES 

When considering densities typical of collapsed objects 
at z = we find that the average density profile of the 
haloes becomes shallower in the core and the core en- 
tropy level increases when a magnetic field is added to 
our simulations. On the other hand, the density profile 
in the outskirts is not affected, even in the presence of a 
strong field because it is mainly in the inner regions that 
the magnetic pressure can reach equipartition with the 
thermal pressure (see Figure 2), while in the outskirts it 
is always negligible. This shallowing of the gas density 
profile results in a more extended gas distribution and 
naturally leads to a reduction in the baryon fraction 
within our objects in the presence of a magnetic field 
(Figure 8) (see Ettori et al. (2006) for a wider discus- 
sion of this topic). For a reasonably high field strength 
and slope (run B0960) the baryon fraction is reduced to 
60% of the cosmic value, with the effect becoming more 
dramatic as the halo mass decreases. 

Figure 9 displays the ratio of the gas mass within 
r 2500 (the radius that encloses an overdensity of 2500) 
for the models with various magnetic field strengths 
compared to the gas fraction obtained in the run with- 
out a magnetic field. For regions with overdensities ex- 
ceeding 2500 times the mean (which corresponds to the 
halo core region which is typically well observed in X- 
rays) the baryonic mass fraction can be even more dra- 
matically reduced, although the effect is only particu- 
larly pronounced for higher field strengths. 
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Figure 8. Gas fraction within the virial radius as a function of 
the virial mass at z = for a run with no magnetic field and with 
B = 10" 9 G and a = 1. 



5 SUMMARY & CONCLUSIONS 

We have presented a simple approximate scheme that 
implements an isotropic magnetic field within standard 
SPH. This method allows various field strength and am- 
plification schemes to be quickly and efficiently tested. 
Using our approach we have modelled the full range of 
observed field strengths and shown that the stronger 
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Figure 9. Relative gas mass within r25oo (the radius enclosing 
an overdensity of 2500) in runs with a magnetic field versus the 
total virial mass in the run without a magnetic field. The different 
symbols mark different runs. 



fields can have a significant affect on the halo core re- 
gions. 

It is important to note that this work could equally 
well be applied to any physical process that supplements 
the usual gas pressure support with additional pressure 
which has a power-law relationship to the density. This 
could include local turbulence (Dolag et al. 2006), the 
injection of hot gas from supernovae (Silich, Tenorio- 
Tagle & Anorve-Zeferino 2005), radio bubbles (Ensslin 
2003; Soker & Pizzolato 2005; Fabian et al. 2006), or the 
activity of galactic nuclei (Sijacki & Springcl 2006). As 
would perhaps have been expected smaller haloes are 
more significantly altered by these processes, particu- 
larly when there is a strong density dependence. 

Any physical process that effects the central bary- 
onic density within dark matter haloes is astrophysi- 
cally interesting because such a process will naturally 
alter star formation rates, feeding of any central black 
hole and the global X-ray emission from the object (to 
name but a few). Even a small lowering of the cen- 
tral density can have a dramatic effect because many 
of these associated processes operate on a high power 
of the local density (for instance, the X-ray emission, 
largely due to thermal bremsstrahlung, is related to the 
density squared). 

If it indeed turns out that large scale astrophysi- 
cal magnetic fields are similar in strength to those sug- 
gested by Dolag et al. (2001) then in the main they 
will have little effect, particularly if they remain un- 
tangled and so amplify with the density to a low power. 
If the magnetic field in collapsed regions becomes tan- 
gled however, and so grows with a steeper power law 
dependence on density then their affects could be sig- 



nificant, particularly within the very central regions of 
small haloes. 
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